HUPD-9403 
March, 1994 



Computation of the Heavy-Light Decay Constant using 
Non-relativistic Lattice QCD 

Shoji Hashimoto 

Department of Physics, Hiroshima University 
Higashi- Hiroshima, Hiroshima 724 
Japan 



abstract 

We report results on a lattice calculation of the heavy-light meson decay 
constant employing the non-relativistic QCD approach for heavy quark and 
Wilson action for light quark. Simulations are carried out at (3 — 6.0 on a 
16 3 x 48 lattice. Signal to noise ratio for the ground state is significantly 
improved compared to simulations in the static approximation, enabling us 
to extract the decay constant reliably. We compute the heavy-light decay 
constant for several values of heavy quark mass and estimate the magnitude 
of the deviation from the heavy mass scaling law fpy/mp = const. For the 
B meson we find fs = 171 ±22^ MeV, while an extrapolation to the static 
limit yields f s B tatic = 297 ± 36l^ MeV. 



1 Introduction 



Computation of electro-weak transition matrix elements of hadrons is one 
of the most important goals of numerical simulations of Lattice QCD. In 
particular the decay constant of heavy-light mesons which consist of one 
heavy quark (Q) and one light anti-quark (q) has great phenomenological 
and theoretical interest in that finding its value represents a crucial step to 
extract the Cabibbo-Kobayashi-Maskawa matrix elements from experiment. 

For heavy-light mesons such as B and D the method of heavy quark 
effective theory [Q] is applicable. In the static limit in which the heavy quark 
field is replaced by its infinite mass limit, implementation of the heavy quark 
effective theory on the lattice is straightforward and was first proposed by 
Eichten0 and by Lepage and Thacker0. 

The initial attempt to simulate heavy-light systems in the static limit 
was made by Boucaud et al.||. They found that the signal of the heavy- 
light meson correlation function is extremely noisy, which made it impossible 
to identify a plateau in the effective mass corresponding to the ground state 
meson. In order to avoid this problem several groups || [7], ^ |[| applied the 
"smeared" source technique, with which the signal for the ground state is 
expected to be enhanced. These groups reported large values of fs in a range 



j static _ 300-500 MeV. However, it was pointed out by two groups |10|, |TT 



that the extracted values of Jb strongly correlate with the size of regions over 
which the source of the heavy-light meson is smeared. This systematic effect 
casts doubt on the reliability of the results from the static approximation. It 
is most likely that one has failed to identify the correct ground state when 
one used arbitrarily chosen smeared sources. 

The reason why large statistical fluctuations occur when one applies the 
static approximation has been explained by Lepage [|l2f with a very simple 
way of estimating the statistical noise of the hadron correlation functions. 
He pointed out at the same time that the use of heavy quark effective theory 
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including the kinetic term should decrease the statistical noise in correlation 
functions. Taking into account the kinetic and other 1/m.Q < terms where 
rriQ is the heavy quark mass is also necessary to estimate the size of 1/tuq 
corrections in physical quantities. This in fact is very important to extract 
physical predictions for the B meson decay constant. 

The heavy quark effective theory including l/rug terms is called non- 
relativistic QCD. In this paper we report results of a lattice calculation of 
the heavy-light decay constant using non-relativistic QCD for heavy quark 
and Wilson action for light quark. Our simulations are carried out in the 
quenched approximation on a 16 3 x 48 lattice at (3 = Q/g 2 = 6.0. We examine 
to what extent the inclusion of 1/m.g terms improve signals for the ground 
state over those in the static limit, and check the independence of results for 
the decay constant on the smearing size. Employing a set of heavy quark 
masses in the range rriQa = 10 — 2.5, we examine the magnitude of deviation 
from the heavy quark scaling law fp^Jrrip =constant. Combining the results 
of this analysis with an estimate of the renormalization factors including 
tadpole improvement ||14|| , we extract physical prediction for the B meson 
decay constant. 

This paper is organized as follows. In sect. 2 we briefly review the for- 
mulation of non-relativistic lattice QCD. In sect. 3 we describe details of 
our simulations and present our results for the signal to noise ratio and the 
heavy- light decay constant. In sect. 4 we discuss perturbative determination 
of renormalization factors, and present our results for the physical value of 
the B meson decay constant. Our conclusions are given in Sec. 5. 
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2 Non-relativistic Lattice QCD 



The general form of the Lagrangian of heavy quark effective theory is given 
as a power series in the inverse heavy quark mass thq 1 : 

£ = £ + A + (1) 
The first two terms are given by 

£„ = QHD Q (2) 

where Q(x) is a two component heavy quark field, the color SU(3) co- 
variant derivative, and B denotes the chromomagnetic field. The leading 
order term £ represents the infinite mass or the static limit. The next to 
leading order term Ci consists of the non-relativistic kinetic energy and the 
Pauli spin-magnetic interaction for heavy quark. For studies of heavy-heavy 
systems (\l/'s and T's) a power counting rule in terms of the heavy quark 
velocity instead of the inverse heavy quark mass is applicable ]T^]. Since the 
spin-magnetic interaction term gives higher order effects in the heavy quark 
velocity, one can omit this term. We are, however, interested in the heavy- 
light system for which the heavy quark velocity is not a good expansion 
parameter. One therefore must properly include the spin-magnetic interac- 
tion term, and we treat the heavy quark effective lagrangian including all 
1/rriQ terms. 

The lattice action we use for simulations is given by 

S Q = £ + # (n) + c^-cr ■ B]Q{x) (4) 



xt 



where 



# (n) = l-(l-^) n , (5) 
3 A A ■ 

H = (6) 
j=i Zm Q 
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The lattice spacing is denoted by a and the gauge covariant difference oper- 
ator is defined by 

A M Q(x) = U li {x)Q(x + A) - Q( x ) 
A^Q(x) = Q(x)-Ul(x-fl)Q(x-fl) (7) 

where U^ix) is the lattice gauge link variable. The action with n > 2 
represents a modification in order to stabilize high frequency modes for 



rriQ < 3/n[|13l- The spin-(chromo)magnetic interaction term is included for 
keeping consistency of the I/tuq expansion even though it is not expected to 
give much effect on the pseudo-scalar decay constant. The usual clover- leaf 
definition is employed for the chromo- magnetic field operator B. For the 
coupling of the spin-magnetic interaction term we use the tree level value 
c=l. 

The 4-component Dirac field of the heavy quark is expressed in terms of 
the 2-component field Q(x) as 

*M = (^ CT q % w )+o ( iM). (8) 

We omit the lower components for heavy quark for simplicity and define 
bilinear operators composed of the heavy and light quarks as 

O r (x) = ( Q^(x), )Vq(x). (9) 
For instance the lattice axial vector current is given by 

A,(x) = ( Q*(x), h, l5 q(x). (10) 
The heavy quark propagator is obtained solving the evolution equation 
G(x,t+1) = U ] xM [(l - iff) n G(a:,t) + $(aj)<y] (11) 

where $(cc) is the source function of the heavy quark. Computation of this 
deterministic evolution equation is much faster than solving the inverse of 
the Wilson quark operator. 
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3 Simulation 



3.1 Simulation Parameters 

Our numerical simulation is carried out with 40 quenched configurations on 
a 16 3 x 48 lattice separated by 1000 sweeps of the 5-hit Cabbibo-Marinari- 
Okawa algorithm fll5[ at (3 = 6.0. 



The heavy quark masses used are mqa = 1000, 10.0, 7.0, 5.0, 4.0 with the 
n = 1 action and 5.0,4.0,3.0,2.5 with n = 2. The values mQ=5.0 and 4.0 
are used for both n = 1 and 2 actions in order to check consistency of results 
for the two actions. The large value tuq = 1000 is taken to compare with 
results of the static approximation. 

For the light quark we used the Wilson action with the hopping parameter 
K = 0.1530, 0.1540 and 0.1550. The critical hopping parameter is K c = 
0.15708(2). The 7r and p meson masses and the pion decay constant we 
obtained using the standard procedure are given in Table |l|. The values of m p 
and f n extrapolated to K — > K c enables us to set the lattice spacing in terms 
of the physical values m p =770 MeV and / 7r =132 MeV, yielding a~ 1 =2.3(3) 



GeV consistently for both m p and f w where we used Z A = 0.86 fig] 



3.2 Signal to Noise Ratio 

For each heavy and light quark masses we have measured the local-local 
correlation function 

C LL (t) = J2(0\TA,(x)Al(0)\0) (12) 

X 

where (x) is the heavy-light axial vector current given in eq. ([T0|) . For large 
enough times this correlation function is dominated by the ground state, i.e. 
the heavy-light meson of lowest mass, 

C LL (t) — { -^^~e- Et (13) 
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where E is the binding energy for the ground state. 

In Fig.^ we plot the effective binding energy defined by 

E ef f(t) = -\o g (^-^), C(t) = C LL (t) (14) 

for m,Q=1000 (open circles) and 5.0 (filled circles). It is clear that the signal 
is far better for mq=5.0 for which we find a clear plateau at t > 12. It is 
well known that the signal for the correlation function in the static limit is 
very noisy, which we also observe here for tuq = 1000. 

The improvement of ground state signals for large but finite values of niq 



can be qualitatively understood from the estimate of the relative error JT2" 
8C(t) 



oc exp 



{E{Qq) _ E(QQ) +m7T)t 



(15) 



C(t) 

where E(Qq) and E(QQ) are the binding energies of heavy-light and heavy- 
heavy mesons respectively. Values of binding energies for K = 0.1530 are 
listed in Table 0. For finite itlq the negative contribution from E(QQ) sig- 
nificantly reduces the value of the exponential slope from that in the static 
limit where E(QQ) = 0, leading to a much milder growth of the relative 
error. 

We found that our data for 5C(t)/C(t) are quantitatively consistent with 
the above estimate. In Fig.^| we show typical examples of the relative error 
5C(t)/C(t) where solid lines indicate the slope expected from the measured 
values of the binding energies and according to (|i~5"l) . Fitting 5C(t)/C(t) 
with the exponential function exp (at) we compare the results for a with 
the estimate from the binding energies in Table 0. We observe quantitative 

oy 
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agreement between a and E(Qq) — g (QQ) +m ' r except for a few cases. 



3.3 Smearing 



In ref.|]T0[ we reported a disappointing fact that the measured values of the 
heavy-light decay constant in the static limit actually depend on the choice of 
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smearing of the axial vector current. Here we study if this problem is avoided 
in non-relativistic QCD. We use the cube smearing as a typical choice of the 
smearing function and examine the dependence of the decay constant on the 
size of the cube. The smeared current is defined as 

A U X ) = 4s E^fc), 0)7^75<?Or) (16) 
n i 

where the sum is over points contained in a cube of a size n 3 centered at x. 
We employ the Coulomb gauge fixing instead of inserting gauge links between 
the heavy and light quark fields. We compute the local-smeared (LS) and 
smeared-smeared (SS) correlation functions defined by 

C LS {t) = £<0|TA 4 L ( M )Af(0,0)|0> (17) 

X 

C ss (t) = J2(0\TAl(x,t)Ai(0,0)\Q) (18) 

X 

for the sizes of smearing 3 3 , 5 3 , 7 3 and 9 3 . Fitting these correlation functions 

as 

C LS (t) — ► Z LS exp{-Et) (19) 
C^ft) > z s/L (20) 

for large enough t regions, we calculate the heavy-light decay constant using 

1 Z LS 

-f P m P = (21) 

where Za is the renormalization factor of the lattice axial vector current. 

In Fig|| we plot the effective binding energies of the LS (filled circles) 
and SS(open circles) correlation functions with the 5 3 smearing. At wq= 5.0 
(Fig.|3](a)) clear signal of the ground state is observed beyond t ~ 10 both 
for LS and SS correlators. Furthermore the values of the binding energies 
are consistent. At mg=1000 (Fig.|](b)), on the other hand, it is impossible 
to identify a plateau. 
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In Fig.[|we plot the values of fp^Jmp/ ' Za extracted from fits of the cor- 
relation functions over 4 time slices t m i n <t< t max = t m i n + 4 for various 
smearing sizes. For each group of data points t min is taken to be t min =6, 
8, 10 and 12 from left to right. At mg=5.0 (Fig.f|(a)) we observe that the 
estimates converge to the same value after t min ~ 10 for all the smearing sizes 
including the case of no smearing (denoted by l 3 in Fig{|(a)). This gives us 
confidence that the asymptotic region where the ground state dominates the 
correlation function is reached at t ~ 10-12. Furthermore it is interesting to 
notice that the magnitude of errors are similar for various smearing sizes in- 
cluding the case of no smearing. This indicates that the smearing technique 
does not improve statistics for this case. 

At mQ=1000 (Fig.f|(b)) the situation is quite different. Here only results 
for t m i n =6 and 8 are available because of rapid growth of noise at t > 10. For 
these fitting intervals the data still depend on the size of smearing, showing 
that the asymptotic region is not yet reached. It is essential for calculations 
in the static limit to use some method which enhance ground state signals 
in the region t < 10. We do not use the data at mq=1000 in the following 
analysis. 

3.4 Decay Constant 

Because the results for the decay constant for rriQ < 10.0 do not depend on 
the smearing size for t m i n > 10, we choose the cube smearing of size 5 3 and 
extract fp^/mp/ Za from a global fit over the interval 10 < t < 20. Other 
choices give similar results. 

In Table |3| we summarize our results for the binding energy E and the 
decay constant fp-\Jmp/ZA at each value of the hopping parameter of light 
quark K and the heavy quark mass niq. Errors are estimated by the single 
elimination jackknife procedure. We observe that the two actions with n = 1 
and 2 yield consistent values for fp-^Jmp/ Za for tuq = 4.0 and 5.0 where 
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both actions are employed. For each tuq we extrapolate the results at three 
values of K linearly in 1/K to 1/K C (see Fig.[5]), the results of which are also 
given in Table [| 



4 Extracting Physical Value of 

In order to estimate the physical value of the B meson decay constant fs we 
have to determine the heavy quark mass mq corresponding to the b quark and 
the renormalization factor of the axial vector current A complete one- 
loop perturbative calculation necessary for this purpose is not yet available 
for our non-relativistic QCD action (HI) with c—1. Davies and Thacker[f2l|, E2 



however, have reported the one-loop results for the case of c=0, and we use 
their results incorporating the tadpole improvement of lattice perturbation 
theory EF 



4.1 Improved Perturbation Theory for Non-relativistic 
QCD 

To one-loop order the inverse heavy quark propagator Aq can be written in 
the form, 

,2 



Ac 



(e ipo - 1 - Ag ) + (1-Bg 



2n P 



2m 



Q 



(22) 



where A, B and C are numerical constants. We define the mass renormal- 
izaion factor Z m , the energy shift Eq and the wave function renormalization 
factor Zq by 

dpo 1 



G pert (p, t) 



Jpot 



n 27T A, 
1 



(23) 



rypert 

Zq exp 



1+Ag 2 )-{1-Bg 



2n P 



2m Q 



t-i 



-t I E p Q ert + 



P" 



2Zm rt niQ j _ 
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We then obtain 



rppert 
ypert 



rypert 



-ln[l + V] 

1 + (A + B)g 2 
1 



l-(C-A)g* 

For these quantities the mean-field improved |23[] expressions are given by 

3 



(24) 



F rnf 



-In 



U Q 1 



nrriQ 



zz f = - i 

MO V 



1 



nrng 



Uq 



(25) 



where «o is a mean-field value of the lattice link variable U^{x). As a def- 
inition of Mo we take the simplest choice uq = (^Trf/p^) 1 / 4 for which the 
perturbative expansion is given by u p ert = 1 — j^g 2 - 

Tadpole-improved one-loop results for the renormalization factors are ob- 
tained by combining (^5|) and (^4j) after removing one-loop terms of the 
mean-field expressions ([ 



r-itad—imp 



hi 



from A, B and C. This gives 
3 



u 1 - 



nm Q 



;i-«o) {i + ^ 2 } 



^tad—imp | j 

ytad—imp 



nm Q 
1 



(l-«o)) 



'l-u )){l + (A + B)g 2 } 



-{l + (C-A)g 2 }. 



(26) 



where A, (A + B), (C — A) are again numerical constants defined by 



.4 



A + — 1 + — 



12 
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(A + B) = (A + B)--j-(l- J> 



12 y nrriQ J 
(C~A) = (C-A)--U (27) 

In Table ^ we summarize numerical values of these quantities for some typical 
values of ttiq where we use the values obtained by Davies and Thacker for the 
numerical coefficients A, B and C . Smallness of the 'tilde' quantities com- 
pared with the original values demonstrates that the tadpole improvement 
works well for these quantities. 

The values of E , Z m and Zq are also tabulated in Table ||. For the 
coupling constant g 2 we take g\r{ji /a) = 1-96 at (3 = 6.0[]23]]. In principle it is 



more desirable to use gy(q*) with a properly determined scale q* as proposed 
in [ 23] . Morningstar has calculated this scale q* for a more complicated non- 
relativistic QCD action|23]] and obtained g*=0.67a~ 1 for Z m and 0.81a -1 
for E . The corresponding values of gy(q*) (it e.g., g v (0.67/a) = 3.41) are 
substantially larger than gy(7r/a). Because of the smallness of the tadpole- 
improved one-loop coefficients, however, the renormalization constants are 
modified only slightly; the combination Z m mQ — Eq which is relevant for 
the B meson mass (see below) changes only by 2-3%, and the change of the 
heavy quark wave function renormalization factor Zq is also at the level of 
a few percent. 

4.2 b Quark Mass 

The heavy-light meson mass is given by 

m P = Z m m Q - E + E Qq . (28) 

In Fig. || we plot our numerical result for the binding energy Eq^ as a function 
of 1/rriQ. The solid line Eq^ = 0.60 + 0.33/mQ fits the data very well. Using 
this fit and the perturbative result for Z m and E discussed in Sec. 4.1, 
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together with the inverse lattice spacing a -1 = 2.3(3) GeV obtained from 
the p meson mass, we find mp = 12.6(1.6), 7.2(9) and 5.4(7) GeV for rriQ 
= 5.0, 2.5 and 1.8 respectively. This shows that the heavy quark mass of 
tuq = 1.8 in lattice units approximately corresponds to the b quark. This 
value is consistent with the estimate given by Davies and Thacker using T 



spectroscopy p5[. In the following analysis of the B meson decay constant 



we use rriQ = 1.8 for the b quark mass. 

4.3 Axial vector current renormalization factor Za 

The renormalization factor Za for the axial vector current can be written as 



Z 



cont 
A 



Za = t&s-. (29) 

L A 

Here Z A ont is the renormalization constant for the continuum axial vector 
current renormalized with the MS scheme [[26]] given by 

ZT 1 = i + Y^9 2 (~ M^r ) 2 - 1) (3°) 

where we use m C Q nt = Z m niQ as the heavy quark mass. We ignore infra-red 
divergent logarithms since it cancels with the lattice counterpart. The lattice 
renormalization factor Z^ tt can be written as 

Z latt = £V2 Z 1/2[1 + ylatt g 2} (3^ 

where Zq is the wave function renormalization factor for the heavy quark 
discussed in Sec. 4.1. For the Wilson light quark our field normalization 
includes the conventional factor \[TK C . The tadpole-improved one-loop result 
for the wave function renormalization factor Z q is then given by|23[ 



2K 

= — | [1 + 0.0043/]. (32) 



Finally the vertex correction V latt is not affected by the tadpole improvement. 



We use the value given in ref. 22 1 for the spinless non-relativistic QCD action 
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(c=0) since theirs is the only result available for this quantity at present. 
Combining these quantities we obtain the final expression, 

Z a = fl-J-(l-^o)V /2 f^V /2 (33) 



x 



nm,Q I \2k c 



where D is a numerical constant. 

In Table f| we list representative values of V latt ,D, 1/Z^ a and Z A - It is 
interesting to note that the value of Z^ tt is almost independent of 1/rriQ while 
Za shows logarithmistic dependence and even diverges at rriQ = oo reflecting 
the form of Z^ nt (eq. (|30|) ) . We also observe that the expansion coefficient 
D is still large after the tadpole improvement, suggesting the possibility 
that higher order coefficients may not be small M. The main contribution to 
the large coefficient is the vertex correction V latt which is unmodified by the 
tadpole improvement. This is a general feature for renormalization constants 
of bilinear operators of light-light, heavy-light and heavy-heavy quarks. 

The large one-loop coefficient increases the magnitude of uncertainties 
coming from the choice of the scale q* in the coupling constant gy{q*). In the 
static limit, the optimal value of q* is estimated |27j as q* = 2.18/ a for which 



g v (q*)=2.22 at (3 — 6.0. This choice of q* leads to Z A =0.QQ1 at m Q = oo, 
while another choice q* = 1/a for which g v (q*)=3.10 gives Za=0.570 which 
is 14 % smaller. In the following we use q* = 2.18/ a for a reference value, 
keeping in mind this large systematic uncertainty. 



4.4 Decay Constant 

Our results for fp^/mp/ZA after extrapolation to K = K c for light quark 
are plotted as a function of I/tuq in Fig.[7] (see Table |3| for numerical values). 
Circles and triangles are for the results obtained with the n = 1 and 2 actions 
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respectively. The two actions yield consistent values for ttlq = 4.0 and 5.0 as 
we already noted in Sec. 3.4. 

We observe in Fig. [7| a clear deviation from the scaling law of the heavy- 
light decay constant in the heavy mass limit given by 

fpyjmp/ Za = const. (34) 

In order to evaluate the magnitude of the deviation we fit the data with the 
form 

f P ^/z A = u P V^p~/z A y taUc (i - — ) . (35) 

V m Qj 

In order to estimate systematic uncertainties in the fit due to a slight curva- 
ture in the 1/mq dependence of fp^/mp/ZA, we list in Table |] the results 
of fit for four representative selection of data points. The first two choices 
employ all data points except at itlq = 4.0 and 5.0 where the results with the 
n = 2 action is chosen for the choice (a) and those with the n = 1 action for 
the choice (b). The choice (c) uses data for 5.0 > itlq > 2.5 with the n = 2 
action and (d) for 10.0 > itlq > 4.0 with n — 1. Errors given in Table |B| are 
estimated by a single elimination jackknife procedure. 

For the decay constant in the static limit ( fp s/Tn p / Z A ) staUc all four choices 
yield values consistent within 10-15%. As a representative value we quote 
the result for the choice (a); 

{fpV^P~/Z A ) staUc = 0.292(35), c = 1.04(44) (36) 



The extrapolated value (pq) is compared with the results of other groups 
obtained at (3 = 6.0 using the static heavy quark propagator in Table 0. The 
results are consistent in view of the systematic uncertainties in the static 
results depending on the detail of smearing. 

Using the value fl3"6"| ) for (fpy/m p/Za) staUc we obtain the heavy-light decay 
constant in the static limit, 

-1 \3/2 



(fpV^Y tattc = 0-682 ± 82tf [^) GeV 3 / 2 (37) 
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where the first error is statistical and the second one is systematic reflecting 
variation of central values among the four selection of data points in Table || 
For the renormalization factor of the axial current Za which diverges in this 
limit we use the value evaluated at the 6-quark mass as a normalization as 
usual. For the B meson decay constant evaluated in the static limit this 
yields 

jstatic = 2Q7 ± 3Q+15 f ^A_\ ( A 1 \ M y /gg) 

JB ~ 30 V0.67/ V2.3GeV/ v ; 

In order to obtain the B meson decay constant including 1/rriQ corrections 
we extrapolate the fit ( p5|) to the b quark mass of rriQ—1.8 as was discussed in 
Sec. 4. 2, employing the single elimination jackknife procedure for estimating 
the statistical error of the extrapolated value. The result shows a sizable 
variation depending on the selection of data points. We quote the value for 
the fitting choice (d) in Table |6] employing data for 5.0 > rriq > 2.5, which 
are close to the extrapolated point uiq = 1.8 and hence should be more 
reliable; 

f B = 171± 22- MeV. (39) 

where the first error is statistical and the second one is systematic showing 
the uncertainty originating from the choice of data points for the fitting. 
It should be remembered that a systematic uncertainty also exists in the 
determination of Za and a -1 . 



The UKQCD collaboration |28[ recently reported a preliminary non-relativistic 
QCD result for obtained in a simulation carried out at uiq = 1.7 on a 
16 3 x 48 lattice. Their result fs yfrriB / Z a = 0.16(3) is consistent with our 
value obtained by an extrapolation in 1/tuq. 

It is also interesting to compare our results with the results obtained 
using propagating quark for the heavy quark. Two groups have reported the 
results of fp-s/mp obtained using the propagating quarks at (3 =6.0. In Fig.|| 
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we plot the quantity 



( 



a s (m P ) 
a s (m B ) 



) 



2/11 



F P 



fp^/mp 



(40) 



as a function of 1/rrip. The normalization factor (a s (mp)/a s (mB)) 2 ' n is 
introduced in order to absorb the logarithmistic divergence of fp^Jmp at 
rap = oo originating from Z^ nt (eq.(^)). As a coupling constant we use 
ay with Ay=0.169 at /5=6.0. Open squares and tringles are for results of 
the PSI-Wuppertal collaboration using the Wilson action with the stan- 
dard normalization \J2K (triangles) and with the improved normalization 



1 — 3K/4K C (squares). Open circles and diamonds are for results of the 



Closed symbols are for results of this work. We can see that our results are 
consistent with the results of the clover action in view of the 1/mp depen- 
dence of Fp. For the Wilson action there is a source of large systematic 
uncertainty in the choice of the normalization for this heavy mass region and 
our results are not seen to be consistent with both of these choices. 

5 Conclusions 

In this article we have reported on a calculation of the heavy-light decay con- 
stant using non-relativistic lattice QCD. We found that ground state signals 
in the correlator is significantly improved compared to those in the static 
limit, and that the degree of improvement of signal to noise ratio is in a 
quantitative agreement with the estimate of Lepage in terms of binding en- 
ergies. As a result we could extract properties of the ground state reliably. 
In particular an apparent dependence of the decay constant on the size of 
smearing for source, which affected previous attempts in the static limit, is 
absent. 

Our result for the B meson decay constant shows that the 1/mg correction 




(clover) fermion action. 
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to the static limit is quite significant even for the b quark. This points 
toward the necessity of a more complete calculation to order 1/mg than was 
attempted here, and eventually a calculation including terms, for a 

precise determination of the B meson decay constant. The improvements of 
the present work needed to order 1/mg are the inclusion of 1/mg terms in 
the axial vector current and one-loop calculation of renormalization factors 
including the spin-magnetic interaction. It is also desirable to estimate two- 
loop vertex corrections in view of the large one- loop coefficient. We leave 
these problems for future investigations. 



Acknowledgements 

I am grateful to A. Ukawa for many discussions and a critical reading of 
the manuscript. I would also like to thank A. Kronfeld, P. Mackenzie and 
M. Okawa for useful discussions. The numerical computation were performed 
on HITAC S820/80 at KEK. I am grateful to the Theory Division of KEK 
for warm hospitality. The author was supported in part by the Grant-in-Aid 
of the Ministry of Education under the contract No. 040011. 



17 



References 

[1] For a recent review, see B. Grinstein, and M. Neubert, SLAC-PUB-6263 
(1993). 

[2] E. Eichten, in " Field Theory on the Lattice", Nucl. Phys. B (Proc. 
Suppl.) 4 (1988) 147. 

[3] G.P Lepage and B. Thacker, ibid, 199. 

[4] For a recent review, see C.W. Bernard, review at lattice 93 to appear in 
the proceedings, hep-lat. 9312086. 

[5] Ph. Boucaud et al., Phys. Lett. B220 (1989) 219. 

[6] C. Alexandrou, F. Jegerlehner, S. Giisken, K. Schilling and R. Sommer, 
Phys. Lett. B256 (1991) 60. 

[7] C.R. Allton, C.T. Sachrajda, V. Lubicz, L. Maiani and G. Martinelli, 
Nucl. Phys. B349 (1991) 598. 

[8] E. Eichten, G. Hockney and H.B. Thacker, in Lattice 90, Nucl. Phys. B 
(Proc. Suppl.) 20 (1990) 500. 

[9] C. Bernard, J. Labrenz and A. Soni, Nucl. Phys. B (Proc. Suppl.) 20 
(1990) 488. 

[10] S. Hashimoto and Y. Saeki, Mod. Phys. Lett. A7 (1992) 387. 

[11] C. Bernard, CM. Heard, J. Labrenz and A. Soni, Nucl. Phys. B (Proc. 
Suppl.) 26 (1992) 384. 

[12] G.P. Lepage, in "Lattice 91", Nucl. Phys. B (Proc. Suppl.) 26 (1991) 
45. 

[13] B.A. Thacker and G.P. Lepage, Phys. Rev. D43 (1991) 196. 

18 



[14] G.P. Lepage et al., Phys. Rev. D46 (1992) 4052. 

[15] N. Cabbibo and E. Marinari, Phys. Lett. B119 (1982) 387; M. Okawa, 
Phys. Rev. Lett. 49 (1982) 353. 

[16] G. Martinelli and Y.C. Zhang, Phys. Lett. B123 (1983) 433. 

[17] APE collaboration, Nucl. Phys. B (Proc. Suppl.) 30 (1993) 469. 

[18] C.W. Bernard, J.N. Labrenz and A. Soni, preprint UW/PT-93-06, hep- 
lat.9306009. 

[19] C. Alexandrou, S. Giisken, F. Jegerlehner, K. Schilling and R. Sommer, 
preprint CERN-TH 6692/92, hep-lat. 9211042. 

[20] UKQCD collaboration, Phys. Rev. D49 (1994) 1594. 

[21] C.T.H. Davies and B.A. Thacker, Phys. Rev. D45 (1992) 915. 

[22] C.T.H. Davies and B.A. Thacker, Phys. Rev. D48 (1993) 1329. 

[23] G.P. Lepage and PB. Mackenzie, Phys. Rev. D48 (1993) 2250. 

[24] C.J. Morningstar, talk at Lattice 93 to appear in the proceedings, hep- 
lat.9311034. 

[25] C.T.H. Davies and B.A. Thacker, Nucl. Phys. B405 (1993) 593. 

[26] See, for example, E. Eichten and B. Hill, Phys. Lett. B234 (1990) 511. 

[27] O.F. Hernandez and B.R. Hill, preprint hep-lat.9401035. 

[28] UKQCD collaboration, C.T.H. Davies, talk at Lattice 93 to appear in 
the proceedings, hep-lat. 9312020. 

[29] C. Alexandrou, S. Giisken, F. Jegerlehner, K. Schilling, G. Siegert and 
R. Sommer, preprint DESY 93-179, hep-lat.9312051. 



19 



K 


m n m p U/Z A 


0.1530 
0.1540 
0.1550 


0.419(4) 0.502(12) 0.105(4) 
0.360(5) 0.461(14) 0.096(5) 
0.291(7) 0.416(21) 0.081(6) 


K c 


0.338(40) 0.067(11) 



Table 1: ir and p meson masses and the pion decay constant f w devided by 
the renormalization constant of the axial current in lattice units. The values 
of m p and J^/Za at the critical hopping parameter are obtained by a linear 
extrapolation in 1/K. 



m Q 


E{Qq) 


E(QQ) 


E(Qq) E ^ 2 +m « 


a 


1000 




0.647(18) 





0.438(18) 


0.466(10) 


10.0 




0.713(24) 


0.510(08) 


0.249(24) 


0.155(51) 


7.0 




0.720(16) 


0.638(12) 


0.192(17) 


0.153(30) 


5.0(n= 


= 1) 


0.740(14) 


0.751(12) 


0.155(15) 


0.149(22) 


5.0(n= 


=2) 


0.738(13) 


0.729(13) 


0.164(15) 


0.147(24) 


4.0(n= 


= 1) 


0.759(14) 


0.827(11) 


0.136(15) 


0.125(22) 


40(n= 


=2) 


0.756(12) 


0.804(09) 


0.145(13) 


0.138(22) 


3.0 




0.786(12) 


0.906(14) 


0.124(14) 


0.116(23) 


2.5 




0.810(15) 


0.977(14) 


0.112(17) 


0.101(23) 



Table 2: The heavy-light and heavy-heavy binding energies E(Qq) and 
E(QQ) at K=0.1530 for the light quark. The combination E(Qq) - g(Q ° 2 )+m " 
determines the behavior of the noise (see Eq.(^) in text), a denotes the 
measured value of the exponential slope. 
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Binding 


energy £ 








m Q 


fT= 0.1530 


0.1540 


0.1550 


7^=0.157 


10.0 

7.0 

5.0 

(n=l) 
(n=2) 

4.0 
(n=l) 
(n=2) 

3.0 

2.5 


0.696(15) 
0.711(14) 


0.681(15) 
0.695(13) 


0.668(17) 
0.679(16) 


0.640(21) 
0.647(18) 


0.734(12) 
0.732(12) 


0.717(14) 
0.714(14) 


0.700(14) 
0.698(14) 


0.666(15) 
0.665(15) 


0.755(13) 
0.750(12) 
0.781(15) 
0.806(13) 


0.737(12) 
0.733(15) 
0.763(14) 
0.788(13) 


0.720(14) 
0.716(13) 
0.746(13) 
0.770(12) 


0.685(12) 
0.682(15) 
0.711(13) 
0.735(15) 


Decay constant fpy/mp/ 'Za 


m Q 


0.1530 


0.1540 


0.1550 


AT C =0.157 


10.0 

7.0 

5.0 

(n=l) 
(n=2) 

4.0 
(n=l) 
(n=2) 

3.0 

2.5 


0.346(19) 
0.328(14) 


0.329(20) 
0.310(15) 


0.314(24) 
0.292(17) 


0.282(32) 
0.257(22) 


0.313(11) 
0.319(12) 


0.294(12) 
0.300(13) 


0.276(13) 
0.282(14) 


0.241(17) 
0.246(18) 


0.302(10) 
0.309(11) 
0.295(09) 
0.287(09) 


0.284(11) 
0.291(11) 
0.278(10) 
0.270(09) 


0.266(12) 
0.272(13) 
0.261(11) 
0.253(10) 


0.232(15) 
0.237(16) 
0.227(14) 
0.220(12) 



Table 3: Binding energy E and decay constant fpyjmp/ ' at if=0.1530, 
0.1540, 0.1550. Values at K c =0. 15708(2) are obtained by a linear extrapola- 
tion in 1/K. 
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m Q 


A 


A 


E Vert 




rptad—imp 


oo 


-0.1684 


-0.0851 


0.401 


0.129 


0.310 


5.0 


-0.2075 


-0.0742 


0.522 


0.202 


0.359 


2.5 


-0.2487 


-0.0654 


0.668 


0.277 


0.414 


1.8 


-0.2809 


-0.0586 


0.800 


0.338 


0.460 


Zfn 


m Q 


A + B 


A + B 


ypert 
m 


Z mf 
m 


^tad—imp 


OO 


0.0686 


-0.0147 


1.13 


1.14 


1.11 


5.0 


0.0434 


-0.0101 


1.09 


1.05 


1.03 


2.5 


0.0707 


0.0374 


1.14 


1.05 


1.13 


1.8 


0.0675 


0.0536 


1.13 


1.02 


1.13 


Zq 


mQ 


C-A 


C-A 


rypert 


2® 


rvtad—imp 


oo 


0.0383 


0.0383 


1.075 


1.000 


1.075 


5.0 


0.0673 


0.0173 


1.132 


1.077 


1.114 


2.5 


0.1003 


0.0003 


1.197 


1.161 


1.162 


1.8 


0.1268 


-0.0121 


1.249 


1.235 


1.206 



Table 4: Renormalization parameters for non-relativistic QCD. The n—1 
action is used at niQ=oo and 5.0 while the n=2 action is used at 2.5 and 1.8. 
As the coupling constant we use gy^n/ a) =1.96 at /3=6.0. 



Z A 


m Q 


ylatt 


D 


1 A 


Z A 


oo 


0.1070 


0.1346 


0.638 


oo 


5.0 


0.0972 


0.1143 


0.654 


0.721 


2.5 


0.0920 


0.1006 


0.655 


0.692 


1.8 


0.0897 


0.0921 


0.650 


0.671 



Table 5: Renormalization factors for the heavy-light axial vector current. As 
a coupling constant we use ffy(2.18/a)=2.22[p7| at /3=6.0. 
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IT) /-i 11 






f ' D . / n? id 7, a 

J B \ I* L B I ^ A 




10.0, 7.0 (n=l) 








(a) 


and 
5.0, 4.0, 3.0, 
2.5 (n=2) 


0.292(35) 


1.04(35) 


0.124(44) 


(b) 


10.0, 7.0, 5.0, 
4.0 (n=l) and 
3.0, 2.5 (n=2) 


0.284(31) 


0.61(21) 


0.187(19) 


(c) 


10.0, 7.0, 5.0, 
4.0 (n=l) 


0.307(44) 


1.03(44) 


0.132(52) 


(d) 


5.0, 4.0, 3.0, 
2.5 (n=2) 


0.262(22) 


0.64(21) 


0.168(22) 



Table 6: Results of fitting with the form (|35"D in the text and the value 
of fpy/mp/ZA at rriQ = 1.8 (last column), (a), (b), (c) and (d) label the 
selection of data points. 



Group 


Lattice 


Smearing 


{fpV^/z A y tauc 


APE|t7j 


18 3 x 32 


cube, n=5 


0.368(14) 






n=7 


0.311(14) 


Bernard et al. IE 


24 a x 39 


cube, n=5 


0.369(18) 






n=7 


0.299(10) 






n=9 


0.241(12) 


PSI-Wuppertal(ll 


12 3 x 36 


exponential, gaussian 


0.327(14) 


UKQCDfgg 


16 3 x 48 


Jacobi 


0.298(10) 


this work 


16 3 x 48 


independent (see text.) 


0.292(35) 



Table 7: Comparison of results for (fp^Jmp / ) statlc obtained by several 
groups at /3=6.0. 
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Figure 1: Effective binding energy of the local- local correlation function for 
niQ = 1000 (open circles) and rriQ = 5 (filled circles), both with K = 0.1530 
for light quark. 
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Figure 2: Relative error of local-local correlation function for uiq = 1000 
(open circles) and uiq = 5 (filled circles), both with K = 0.1530 for light 
quark. Solid lines indicate the slope expected from measured values of the 
binding energies and m n according to (]T5|) . 
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1.2 
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Eeff(t) 
1.2 



0.8 







SHI 

r 



^0 



10 



20 



t 



(b) 



Figure 3: Effective binding energy of the local-smeared (LS) and smeared- 
smeared(SS) correlation functions for (a) mg=5.0 and (b) mQ=1000, both 
with i^=0.1530 for light quark. Smearing size is 5 3 . 
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Figure 4: Dependence of fp^/mp/ZA on the fitting range t m i n < t < t min + 4 
for various smearing sizes at tuq =1000 (a) and uiq =5.0 (b) with fT=0.1530 
for light quark. 
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Figure 5: fpy/mp/ 'Za as a function of 1/2K for mQ=10.0 (circles), 5.0 
(squares) and 2.5 (diamonds). The solid line indicates the extrapolation 
to the chiral limit 1/2K C . 



30 



E Qq 

0.8 



0.7 



0.6 




0.2 0.4 

l/m Q 



Figure 6: The binding energy E Q g of the heavy-light meson as a function of 
I/wig Circles and triangles are for the results obtained with the n—1 and n=2 
actions respectively. The solid line represents the fit Eq^ = 0.60 + 0.33/mg. 
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Figure 7: 1/mg dependence of fp^/mp/ Za after the extrapolation to the 
limit K = K c for light quark. Circles and triangles are for the results obtained 
with the n—1 and n=2 actions respectively. 
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Figure 8: Fp = {a s {mp) / a s {mB)) 2 ^ 11 f py/mp as a function of 1/mp after the 
extrapolation to the limit K = K c for light quark. Open squares and tringles 
are for results of the PSI-Wuppertal collaboration [^9] using the Wilson action 
with the standard normalization \/2K (triangles) and with the improved 
normalization ^1 — 3K/ 4K C (squares). Open circles and diamonds are for 
results of the UKQCD collaboration [|20] using the 0(a)-improved (clover) 
fermion action. Closed symbols are for results of this work. 
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